A Nonequilibrium Field Theory Description of the Bose-Einstein Condensate 
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We study the detailed out of equilibrium time evolution 
of a homogeneous Bose-Einstein condensate. We consider a 
nonrelativistic quantum theory for a self-interacting complex 
scalar field, immersed in a thermal bath, as an effective mi- 
croscopic model for the description of the Bose-Einstein con- 
densate. This approach yields the following main results: (i) 
the interaction between fluctuations proves to be crucial in 
the mechanism of instability generation; (ii) there are es- 
sentially two regimes in the fc-space, with a crossover for 
k 2 /2m ~ 2\\ipo\ 2 , where, in our notation, A is the coupling 
constant and |<£o| 2 is the condensate density; (iii) a set of 
coupled equations that determines completely the nonequi- 
librium dynamics of the condensate density as a function of 
the temperature and of the total density of the gas. PACS 
number(s): 03.75.Fi, 05.30.Jp, ll.10.Wx 

The experimental verification of the phenomenon of 
Bose-Einstein condensation in weakly interacting gases 
has boosted a large number of theoretical investigations 
on the dynamics of weakly-interacting dilute gas sys- 
tems [for a recent review, see e.g. Ref. [|l| and refer- 
ences therein]. Current experiments and planned ones 
make it possible to probe different aspects of the Bose- 
Einstein condensate formation, with great control over 
interactions, trapping potentials, etc. Nevertheless, a 
basic problem not yet fully understood is the following: 
given an initial state, how will the condensate evolve with 
time? In special, the time scales for the condensate for- 
mation and its final size are important quantities involved 
in recent experiments with dilute atomic gases Q . 

On the theoretical side, however, only restrict progress 
has been achieved concerning the problems above. Previ- 
ous studies by Stoof || were able to give a qualitative idea 
of the various time scales involved during the condensate 
formation. In fact, they were the first attempts to an- 
alyze the problem from a microscopic point of view, by 
using the Schwinger-Keldysh closed time-path formalism 
(for reviews, see for instance Refs. ) in the quantum 
field theory description of Bose-Einstein condensation. 
Regarding the condensate growth problem, Gardiner et 
al. |(| have used a quantum kinetic theory to construct 
a master equation for a density operator describing the 
state of the condensate, which is equivalent to a Boltz- 
mann equation describing a quasi-equilibrium growth of 
the condensate. 

In this work we will study the quantum field time evo- 
lution of an interacting homogeneous condensate. Al- 
though non-homogeneity is inherent to current experi- 



ments on Bose-Einstein condensation of atomic gases in 
trapping potentials, we believe that a full understanding 
of the time evolution of even the simpler case of a ho- 
mogeneous gas is still lacking. Besides, as pointed out 
by Stoof in ||, the simplest formulations based on ki- 
netic theory do not allow for the observation of a macro- 
scopic occupation of the one-particle ground state, and 
the question of the instability of the Bose gas system in 
the homogeneous case is a nontrivial one. This makes 
its study an interesting problem, which may shed some 
light on the analysis of the systems under experimental 
investigation. 

We consider the simplest model for a nonrelativistic 
complex Bose field, with a hard core interaction poten- 
tial, whose Lagrangian density is given by (throughout 
this work we use units such that h = 1) 



C 



dt 2m 



+ M 0*0-A(0» 2 , (1) 



where the complex scalar field </>(x, t) represents com- 
plex spinless bosons of mass m, and A is the coupling 
constant, related to the s-wave scattering length a by: 
A = Ana/m. In (jl]) we have also explicitly introduced a 
chemical potential fi that produces a constant total den- 
sity of particles (<t>*4>) = n. We also assume that the 
system is coupled to a heat bath environment with in- 
verse temperature j3 = 1/fcsT. 

We may now perform the standard decomposition Q] 
of the fields {4>,<f>*) into a condensate (uniform) part 
(ipo, v?q) and a fluctuation (nonuniform) part (ip, tp*) that 
describes the atoms outside the condensate, as </>(x, t) = 
ip (t) + ip(x,t) and <f>*(x,t) = + <p*(x,i), where we 
have assumed a homogeneous condensate. Note that we 
take at first tfo(t) as an arbitrary function of time that 
will be determined by the dynamics of the system. 

Substituting the fields above in ([!]), we can readily 
obtain the Bogoliubov spectrum for quasi-particles [§j. 
In particular, the quadratic part of the Lagrangian for 
the fluctuation fields is the Bogoliubov approximation 
for quasiparticles. However, this approximation has a 
flaw: there is no interaction between the fluctuations, 
which is only reasonable at temperatures well below the 
critical temperature for the condensate formation. Our 
aim, in this work, is the study of the condensate evolu- 
tion and, therefore, we must go beyond the Bogoliubov 
approximation. The simplest extension is to implement 
a mean field approximation in the interactions between 
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the fluctuation field. In this way, as we show below, one 
can make clear the appearance of instability modes to- 
wards the condensation formation once the interactions 
between fluctuations are taken into account. With the 
decomposition above, the interaction term for fluctua- 
tions in the Lagrangian becomes X(ip*<p) 2 . The mean- 
field approximation amounts to the following: 

\(tp*<p) 2 =4X(<p*<p)<p*tp + [A(^» 2 - 4\(<p*<p)(p*<p] , 

(2) 

where the first term in the rhs is taken as part of the 
quadratic Lagrangian for fluctuations, and the term in- 
side the square brackets is taken as part of the interaction 
Lagrangian. 

From the decomposition of the fields and (j^), the 
quadratic part of the Lagrangian density for the fluc- 
tuations, C${<p, f*), may be written as 



.d 1 , 

! 1 

dt 2m 



ip + tp* {-\(pl)(p* 



+ <P(-^o )<P 



(3) 



Here, we have used the fact that, under the field de- 
composition in the condensate and out of the conden- 
sate modes, the density constraint then becomes (4>*4>) = 
|</?o| 2 + (y*V) = n - Additionally, assuming that at the 
initial time the system is mostly composed of particles 
outside the condensate, (<f>*4>) — (<fi*f) (at t = 0), sim- 
ple relations involving the generating functional for the 
correlation functions (see, for instance, the last section 
of chapter 2 in 0]) allow us to write the total number 
density n of particles in terms of the chemical potential 
//, valid in the mean-field approximation for the poten- 
tial, as: /it = AXn. Note that this is just the expression 
obtained also in the Hartree-Popov approximation [(Tot , 
which turns out to satisfy the Hugenholtz-Pines relation 
|PH that would be obtained in the equilibrium problem. 
These considerations lead to the quadratic Lagrangian 
for the fluctuations shown above. 

The scenario we have in mind is that for time t < the 
initial state is in equilibrium at a temperature Tj ^> T c . 
At t = the system is then abruptly quenched to a 
much lower temperature Tf < T c . Tf is the temperature 
of the thermal bath in which the system is immersed 
and, of course, it will be the equilibrium temperature 
which the system will reach asymptotically. This kind 
of quench is easily attained in the experiments of Bose- 
Einstein condensation of atomic gases, where the typical 
relaxation time scales are long enough (around ~ 0.1s, 
depending on the temperature [O) to allow for a fast 
drop in the temperature of the system that evolves af- 
terwards out of equilibrium. With this choice of initial 
state, it is reasonable to approximate the dynamics of 
the build up of the condensate state, which at the initial 
time is n con( j.(t = 0) = \<fo(t = 0)| 2 w 0, and the de- 
pletion of the excited states (which at t = it is given 



by n w n cxc .(i = 0) = (ip*ip)) as essentially a two- level 
problem. It is clear that this approximation breaks down 
for temperatures close to the critical temperature, where 
the detailed treatment would require a thorough study of 
the dynamics among the many levels of excited states. In 
the above approximation, the condensate builds up sub- 
ject to the density constraint relation, which may be ex- 
pressed in terms of the averages of the real and the imag- 
inary parts of tp and tp* (ip — <p>\+itp2 and ip* = ip\ — i<p2, 
respectively). Spatial translational invariance yields: 



M*)l 2 

^exc. (t) 



n cxc .(t) = n . 
(pi(%l(t)) 



(4) 



The field averages above can be expressed in terms of the 
Green's functions for <p\ and ip2 as (j — 1, 2) 



/|j^HG>(k,M)] 



(•5) 



where G^ (k, t,t) is defined from the Green's functions 
for the fields in the closed-time path d|| (in momentum 
space) 

G++(k, t, O - G> (k, t, t')9(t - t') + G< (k, t, t')6(t' - t) , 
Gjr (k, t, t') = G> (k, t, t')9(t' -t) + G< (k, t, t')6(t - t') , 
G+-(k,M') = -G< (k,M') , 

Gj j h (k,t,t') = -Gf j (k,t,t') . (6) 

The functions G > and G < satisfy the property 
G^ (k, t, t') — G> (k, t—i(3, t'), which is recognized as the 
periodicity condition in imaginary time (Kubo-Martin- 
Schwinger (KMS) condition). (3 here is the inverse of the 
temperature of the thermal bath and appears here as a 
consequence of the boundary conditions arising from the 
construction of the complex time path. G > and G < are 
constructed from the homogeneous solutions to the op- 
erator of quadratic fluctuation modes, which, using Eq. 
(^) expressed in terms of p\ and p2, are given by (in 
momentum space) 

^ + (£ + »W«)»(M>-o 

The boundary conditions for the solutions of the equa- 
tions above are such that, for t < 0, |<^o(i)| 2 = 0, 
Xi(k, t) = cos(ek^) and % 2 (k, t) = — sin(ekt), where 
£k = k 2 /(2m). In terms of these fluctuations modes, 
the Green's functions are expressed as 

G>(k,M') 



2(1 - e-0 £ *) 

x [xjMxlM+e-^XjMXiM] (8) 
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andG<(k,t,i')=G>(k,i',t). 

By decoupling the set of equations in (0), we can read- 
ily identify that those modes with (k 2 /2m) < 2A|</?o| 2 
are unstable and drive the excited particles towards con- 
densation. Note also that not all the excited particles 
condense, since there will always be a fraction (which 
depends on various parameters for a particular system 
and on the temperature of the thermal bath) of excited 
modes, with high enough frequency, that remains stable. 
This will be clear from our numerical results shown later. 

Using Eqs. ([?]), (||) and the boundary conditions on 
the Green's functions, Eq. (^), together with the initial 
condition on the density (at t — 0, as defined before), 
one can then show that n cxc .(0 can be expressed as (by 
subtracting the zero point divergent contribution) 



.(*) = 



3/2 



d 3 k 



[| Xl (M)| 2 + | X2 (k,i)| 2 ] n k (/3), 

(9) 



where nt(/3) = (e ,3ek — and f3 c is the inverse of 

the equilibrium critical temperature, defined in terms of 
the total gas density, n jl3| It should also be noted that 
in our out of equilibrium approach there are no infrared 
divergences since the finite time is a natural regulator. 
However for the equilibrium t — > oo, the critical temper- 
ature will be modified by the interactions as pointed out 
in Ref. ||. 

The expression above can also be obtained directly 
in terms of the Green's functions for the complex fields 

if, if*: 



<¥>(*)¥>*(*)} = 



d 3 k 



" [-iG^CM,*)] . (ii) 



[-iG> v ,(k,t,t) 



(10) 



In terms of ( JlCj ) and (pd|), we have {<Pi{t)ipi(t)) + 
{<p 2 (t)Mt)\ = «¥>(*)¥>*(*)) + (<p*(tMt)))/2, and the 
KMS condition can be expressed, in this CM SC. clS 
G>,.(M - 0,t') = [G> v (k, *,*')]*, or G> v (k,t- 
0,t') = {G>(k,t,t')}*. 

Equation (ft) is the first order term in the finite temper- 
ature quantum many-body perturbation expansion for 
(if 2 ) . Higher-order corrections for the equal-time two- 
point field averages can be expressed in terms of the coin- 
cidence limit of the (causal) two-point Green's functions 
G vv * and G v * Vl which satisfy the Dyson equations (the 
indices stand for ip and ip*): 



Gij — G°- + G ik Y,kiGi 



3 i 



(12) 



where Ejj is the (matrix) self-energy, and G^- is the 
zeroth-order non-interacting Green's function, satisfying 
the equation (in momentum space) 



dt 



(13) 



One of the advantages of expressing the Green's functions 
in terms of the solutions of (Q) is the possibility of obtain- 
ing, in an unambiguous way, all higher-order corrections 
to the two-point and many-point functions JTq] . 

By using Eq. @ , we can rewrite the constraint on the 
density as 



IvaoC*)! 2 = ~ l I dkk 2 x 



2tt 2 

x [l-(| Xl (k,t)| 2 + | X2 (k,t)| 2 )]n k (/3), (14) 

that reproduces the result obtained by Stoof || for the 
limit t — > oo. Equations (0) and (|l4|) form an integro- 
differential system that may be solved for <fo(t) numer- 
ically, given the initial conditions for <fo(t), X i(fc,i) and 
X2(k,t) mentioned before. Indeed, this system of equa- 
tions determines completely the time evolution of the 
condensate density as a function of the temperature and 
of the total density of the gas. Explicit results for differ- 
ent temperatures are shown in Fig. 1. It is important to 
point out at this stage that the evolution of the conden- 
sate is completely driven by the interactions between the 
microscopic fluctuations of the field around the conden- 
sate. 

Throughout this letter, we have developed an out-of- 
cquilibrium non-perturbative quantum field theory de- 
scription of the condensation process of an interacting 
homogeneous Bose-Einstein gas quenched below the crit- 
ical temperature. In summary, this approach yielded the 
following main results: (i) The interaction between fluc- 
tuations proved to be crucial in the mechanism of insta- 
bility generation; without it, there is simply no macro- 
scopic condensate at all. (ii) There are essentially two 
regimes in the fc-space: for (fc 2 /2m) << 2A|</Jq|; we 
have unstable modes that decay exponentially, while for 
(k 2 /2m) >> 2A|</j 2 |, we have stable modes that oscillate, 
with a crossover for (k 2 /2m) ~ 2A|</? 2 |. (iii) Equations 
(0) and (^) come from a microscopic model for the weak- 
interaction gas, and determine completely the dynamics 
of the condensate. In fact, they are non-perturbative 
and certainly implement a resummation of the ladder 
Feynman diagrams mentioned by Stoof in p). Indeed, 
the highly nonequilibrium character of this description 
should complement the usual approach via Boltzmann 
equation. 

Although we focused this letter on the instability pro- 
cess that generates the condensate (i.e., the short time 
behavior), for t — > oo our results confirm the behavior 
predicted in Ref. || for this limit. However, the equilib- 
rium (t — > oo) values of the condensate fraction are lower 
than the experimental results p2| and the calculations of 
Dalfovo et. al. 0. This may be due to our approximation 
of neglecting incoherent collisional processes, which is a 
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valid approximation in an infinite homogeneous gas at 
very low temperatures and densities but otherwise may 
give an important contribution. We expect that the self- 
consistent inclusion of pair terms should account for most 
of these contributions. We will report on these improve- 
ments of our mean-field approximation in a future pub- 
lication ]l5| ]. 

In spite of the absence of non-homogeneity effects, we 
hope that the approach developed here may be useful 
in the analysis of transients in realistic Bose-Einstein 
condensation experiments with atomic gases. Moreover, 
with a suitable generalization of the formalism presented 
above, we could be able to develop a theoretical descrip- 
tion of the dynamical aspects of a recently proposed ex- 
periment Jlj| (currently in progress) , regarding the Bose- 
Einstein condensation in a weakly-interacting photon gas 
in a nonlinear Fabry-Perot cavity |17|. 
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FIG. 1. Condensate density as a function of time for 
na 3 = 0.01 and Ti/T c = 0.06, T 2 /T c = 0.08 and T 3 /T c = 0.1. 
Here, r = (h/ma 2 )t is a dimensionless time and po = a 3 |c^o | 2 
is a dimensionless density. 
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